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Abstract 

The study of a particle with position-dependent effective mass (pdem), within a double 
heterojunction is extended into the complex domain — when the region within the hetero- 
junctions is described by a non Hermitian VT symmetric potential. After obtaining the exact 
analytical solutions, the reflection and transmission coefficients are calculated, and plotted as 
a function of the energy. It is observed that at least two of the characteristic features of non 
Hermitian VT symmetric systems — viz., left / right asymmetry and anomalous behaviour 
at spectral singularity, are preserved even in the presence of pdem. The possibility of charge 
conservation is also discussed. 
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1 Introduction 



Ever since the pioneering work of Bender et. al. more than a decade ago, the fact that a class 
of non-Hermitian Hamiltonians admit real and discrete spectrum under certain conditions, is 
well established p!| |2] . Non-Hermitian hamiltonians having VT symmetry [P — ?> parity, T — >■ 
time reversal) form a special class in this category, as they admit real and discrete spectrum 
for exact VT symmetry and complex conjugate pairs of energy when this space-time symme- 
try is spontaneously broken, the transition occurring at the so-called exceptional point [3l 0]. 
Naturally, numerous attempts have been made by various scientists to extend the framework of 
quantum mechanics into the complex domain |5i [G] I7j . Theoretical predictions for VT symmetric 
systems exist in quantum field theory, mathematical, atomic and solid state physics, classical 
optics, etc |8j. A pair of coupled active LCR circuits — one with amplification the other with 
equivalent attenuation, exhibits VT symmetry [9]. Of late, experimental confirmation of such 
non-Hermitian VT symmetric concepts have been observed in optics, in VT symmetric crystals 
with a complex refractive index distribution n{x) = no + nji{x) + inj{x), where no represents 
a constant background index, ur^x) is the real index profile (even) of the structure, and n/(x) 
stands for the gain or loss component (odd) [10-17]. Unlike ordinary crystals, complex crystals 
show unique properties — e.g., violation of Fresnel's Law of Bragg scattering, double refraction, 
power oscillations, non-reciprocal diffraction, handedness or left-right asymmetry, anomalous 
transport, unidirectional invisibility, etc. In fact, these unique features together with the occur- 
rence of exceptional points in the discrete spectrum and spectral singularities in the continuous 
spectrum, are characteristic of VT symmetric non Hermitian Hamiltonians, unknown to Her- 
mitian ones. A lasing medium embedded in the spatial region \z\ < ao where the dielectric 
constant satisfies the VT symmetry condition e(— r) = e*(r), behaves as a laser oscillator (LO) 
for positive Im e(f) signifying gain, or as a coherent perfect absorber (CPA) for negative Im e(r) 
representing loss, provided the dielectric constant is real valued and constant e(r) = eo (say) 
outside i.e., for \z\ > qq. While a LO can emit outgoing coherent waves, a CPA can fully absorb 
incoming coherent waves |18j . The interesting part of non-Hermitian Hamiltonians is that even 
a single VT cell can exhibit unconventional features [11] . 

On the other hand the study of quantum mechanical systems with position dependent effective 
mass (pdem) has received a boost in recent times with major developments in nanofabrication 
techniques of semiconductor devices [19-26]. The spatial dependence on the effective mass of 
the particle arises due to its interaction with an ensemble of particles within the device, as the 
particle propagates from left to right. For example, in Al^^Gai-ajAs, as the mole fraction x varies 
along the z-axis, so does the effective mass of the charge carrier (electron or hole). Pdem for- 
malism is extremely important in describing the electronic and transport properties of quantum 
wells and quantum dots, impurities in crystals. He-clusters, quantum liquids, semiconductor 
heterostructures, etc. In a recent work, we obtained the exact analytical scattering solutions 
of a particle (electron or hole) in a semiconductor double heteroj unction — potential well / 
barrier — where the eff'ective mass of the particle varies with position inside the heteroj unctions 
|27j . It was observed that the spatial dependence on mass within the well / barrier introduces 
a nonlinear component in the plane wave solutions of the continuum states. Additionally, the 
transmission coefficient increases with increasing energy, finally approaching unity, whereas the 
reflection coefficient follows the reverse trend, going to zero. 

This study is presented as a sequel to the work done in ref. [27] , in an attempt to extend the 
pdem formalism further into the complex domain. Instances of such attempts are found in other 
works as well — e.g., the problem of relativistic fermions subject to a VT symmetric potential 
in the presence of position-dependent mass was studied in ref. [28], while exact solutions of 
Schrodinger equation for VT / non VT symmetric and non-Hermitian Morse and Poschl- Teller 
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potentials were obtained with pdem by applying a point canonical transformation method in ref 
[29j . In the present work, the main stress will be given on the scattering phenomenon in a VT 
symmetric double heteroj unction with pdem — a special form of semiconductor device consisting 
of a thin layer of VT symmetric material sandwiched between two normal semiconductors, such 
that the mass of the charge carrier (electron or hole) varies with the doping concentration (and 
hence position) within the heterojunctions, but is constant outside. In particular, we shall see 
what new properties (if any) can be expected from VT symmetric heterojunctions with pdem, 
with special emphasis on the behaviour of the reflection and transmission coefficients 

(i) with respect to the direction of incidence of the particle 

(ii) at the spectral singularity. 

Additionally, we shall also derive explicit relations for current and charge densities, to see 
whether these are conserved in such a device. 

For this purpose, we shall consider a double heteroj unction, with the potential function in the 
intermediate region satisfying the VT symmetry condition V{—z) = V*{z), but assuming a real 
valued constant outside 

{Vr{z) + iVi{z) , ai< z < a2 
Vol = VR{ai) , -oo < z < ai (1) 
Vo2 = VR{a2) , 02 < 2; < oo 

where ai,a2 represent the heterojunctions. The mass of the charge carrier is assumed to be of 
the form 

'm{z) , ai < z < a2 

m = I mi = m(ai) , —00 < z < ai (2) 
1712 = m{a2) , a2 < z < 00 

Thus, the mass m{z) and the real part of the potential function viz. Vr{z) are considered to 
be continuous throughout the semiconductor device. We shall mainly concentrate on obtaining 
the exact analytical solutions of the scattering states of a particle with pdem inside a VT 
symmetric double heteroj unction, which is essential to study the nature of the reflection and 
the transmission coefficients. 

The article is organized as follows : For the sake of completeness, the position-dependent-mass 
Schrodinger equation is introduced in Section 2, and the method of obtaining the solutions is 
discussed briefly. To give a better insight into the physical nature of the problem, we shall study 
an explicit model in Section 3, and plot the potential and mass functions as a function of z in 
Fig. 1, and the scattering solutions in Fig. 2. The transmission and reflection coefficients are 
also calculated, and their behaviour is discussed with respect to the relative strengths of the 
coupling parameters of the potential (both real and imaginary parts) and the mass functions. 
Since non-Hermitian Hamiltonians (with constant mass) are known to show peculiar behaviour 
at spectral singularities and also exhibit left-right asymmetry, the transmission and reflection 
coefficients are plotted in Figures 3, 4, 5 and 6, to check whether similar phenomena are observed 
in the presence of pdem as well. Section 4 is devoted to the conservation of charge and current 
densities for non Hermitian VT systems with pdem. Finally, Section 5 is kept for Conclusions 
and Discussions. 



2 Theory 

We start with the basic assumption that the one dimensional time independent Schrodinger 
equation associated with a particle endowed with pdem is the same for Hermitian and non 
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Hermitian systems, and is given by 

Hem{z)^{z) = [TEM{z) + V{z)]xl,{z)=Ei^{z) 

(3) 

V{z) = VR{z)+iVi{z) 
in the intermediate region within the heteroj unctions. Tem is the kinetic energy term given by 

Tem = ^ (jn°'pm^pm'^ + rrCpm^pnf 

1 / 1 



1 V m , 

where p = —ih-^ is the momentum operator. It is to be noted here that the kinetic energy term 

az 

is considered to be Hermitian. The non Hermiticity is introduced through the potential term 
V{z), with an even real part Vr{z) and an odd imaginary part Vj(z). The ambiguity parameters 
a , /3 , 7 obey the von Roos constraint [22] 

a + /3 + 7 = -l (5) 

In the absence of a unique or universal choice for the ambiguity parameters, several suggestions 
exist in literature [301 [311 [32l [33] , 6tc. However, for continuity conditions at the abrupt interfaces, 
and well behaved ground state energy [341 [35] , we shall restrict ourselves to the BenDaniel-Duke 
choice, viz., a = 7 = , (3 = —1. Incidentally, this particular choice consistently produces 
the best fit to experimental results [36]. Furthermore, we shall work in units h = c = 1, and 
use prime to denote differentiation w.r.t. z. Thus, inside the potential well ai < z < 02, the 
Hamiltonian for the particle with pdem reduces to [37] 

whereas, outside the well, z < ai and z > 02, the particle obeys the conventional Schrodinger 
equation : 

^ +Voi,o2] i^{z) = Ei;{z) (7) 



2mi^2 dz"^ 

having plane wave solutions. In case we consider a wave incident from left, the solutions in the 
two regions are 

^PL{z) = e'''^' + Re-''''' , -00 < z < ai 

(8) 

ipR{z) = re*^22 ^ a2<z<oo 
where R and T denote the reflection and transmission amplitudes, and 



fci,2 = ^2mi,2 {E - ^01,02) (9) 
To find the solution in the region ai < z < 02, we make use of the following transformations [38j 
^in{z) = {2m(z)}i/^ ^{p) , p = [ ^2m{z)dz (10) 
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which reduce the Schrodinger equation for pdem, to one for constant mass, viz., 



with 



dp 



V{p) = V{z) 



7 m 



,2 ^, 



(11) 

(12) 



32 rrfi %rn? 

Some definite practical forms of Vi^z) and m(z) give exact analytical solutions of (jllj) . An 
explicit example in the next section illustrates our purpose. 

3 Explicit model : VT symmetric potential well with position 
dependent mass 

We consider the following ansatz for the potential, whose real part describes a diffused quantum 
well 

+ ^ 1 \z\ < ao 



Viz) = { 
Let the mass of the particle be 

m{z) = < 



1 + ^2 1 + ^2 



(13) 



l + «o 



Vr 







\z\ > Oo 



(3' 



2(l + z2) 
2(l + a2) 



, \z\ < ao 



mo , \z\ > uq 



(14) 



where pi, p2, /3 are some constant parameters. 

For a better understanding of the mass dependence and the potential in the semiconductor 
device, we plot m{z) and the real and imaginary parts of V{z) as a function of z in Fig. 1, for 
a suitable set of parameter values, viz., /3 = 4, p = 3, ao = 4. 
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Figure 1: Colour online : Plot showing m(z) and V{z) w.r.t. 



For the spatial mass dependence given by eq. ()14p . eq. (jlOp transforms the coordinate z to 

/9 = /3 sinh-i z (15) 



so that after some straightforward algebra V{p) in eq. (^T2\i reduces to 



V{p) = ^ - ^sech^^ + ipsech ^ tanh ^ (16) 



Thus equation (jlip may be written as 



-r^ + + Visech^/9 — zV2sech p tanhoj (j) = (17) 
dp'' ^ ' 

where = S/?^ - i , p = ^ (18) 

and the parameters Vi, V2 depend on the constants iii.,p,2 and through the equations 

Vi=piP^ + -^ , y2 = //2/3' (19) 
Let us introduce a new variable 

. = (20) 

and write the solutions of ()17p as 

<A = /(l-y)''n(y) (21) 
In terms of the new variable y, equation (jlTh reduces to the hyper geometric equation 



.(1-.)^ + |(^2, + -j-(2p + 2, + l),|- 
- {k^ - (p + g)^}n = 



(22) 



where 



(23) 



Now, (j22p has complete solution |39j 

u = P 2i^i(a,6,c;y) + Qyi-= ai^i (1 + a - c, 1 + 5 - c, 2 - c; y) (24) 

where P and Q are constants, and the parameters a and b are as defined below : 

a = p + q — in , h = p + q + in (25) 

After some straightforward algebra, the final solution to the pdem Schrodinger equation ([3]), 
within the potential well \z\ < uq, is obtained as 

(26) 

+ Q{^^] 2Fi{l + a-c,l + b-c,2-c;y)} 



where y = {l + iz)/2. Outside the well {\z\ > gq), the solutions are given by eq ([8]), with ki = k2. 
The solution in the entire region is plotted in Fig. 2 for the same set of parameter values as in 
Fig. 1, viz., /? = 4, /i = 3, ao = 4. With the help of Mathematica, the constants P, Q and the 
reflection and transmission amplitudes \R\ and \T\ respectively, are determined using modified 
boundary conditions for pdem systems 

1 dip(z) 



the functions il^{z) and 



should be continuous at each heteroj unction iao- That the 



m{z) dz 

effect of pdem is to introduce a non-linear component in the solution is obvious from the figure. 
This finding is analogous to the Hermitian model with pdem |27] . 




Figure 2: Colour online : A plot of Re iIj{z) vs z; Dashed (black) lines show the abrupt heterojunctions 



For real and discrete spectrum, IV2I <Vi + 1/4, implying 



ml < w + ^ 



(27) 



However, in this work we are interested in the scattering states only, i.e., positive . Hence 

1 



E > 



4/32 



Energy 



Figure 3: Plot of |Tp vs E — same for left and right incidence 



The transmission and reflection amplitudes T and R respectively, are plotted in Figures 3, 4 and 
5, for both left and right incidence, with the help of Mathematica. While \T\ comes out to be 
the same for either case, the reflection coefficient depends on whether the particle enters from 
left or right. \R\ is normal {\R\ < 1) when the particle enters from left — the absorptive side 
(Im V{z) < 0), and anomalous {\R\ > 1) when the particle enters from right — the emissive side 
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(Im V{z) > 0). Thus, this phenomenon of left-right asymmetry, characteristic of non hermitian 
VT symmetric potentials with constant mass particles |4H H2t H3] , remains unaltered even when 
the particle mass has a spatial dependence. 




Figure 4: Plot of \R\^ vs E for left incidence 



Energy 



Figure 5: Plot of \R\^ vs E for right incidence 



Another interesting feature worth discussing here is when the well is replaced by a barrier 
i.e., Hi < 0, so that 



V{p) = - + — - ) sech^ p + i/i2/3^sech p tanhp 



(28) 



For particles with constant mass, the potential in ()28p is known to admit a spectral singularity, 
with the reflection and transmission coefficients blowing up at the positive energy [441 145] 



1 



Vo 



1 



when 



V2\>\Vi\ + 



sign of Vi 



(29) 



(30) 



For the pdem non Hermitian heterojunction considered here, the condition (j30p comes out to be 



P2 > Pi + 



1 

2^ 



so that spectral singularity occurs at 



\ip.-p,)^'-l 



(31) 



(32) 



Interestingly, even in a VT symmetric double heterojunction, with a spatially varying mass, 
both and |T|^ blow up at the value of Eg given in eq. ()32p . as observed in Fig. 6. 
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\T\' 



Figure 6: Colour online 
spectral singularity 



Plot of and |Tp vs Energy for VT symmetric barrier : Es shows the 



4 Current and Charge Conservation 

Since we are dealing with non Hermitian VT symmetric systems, it would be interesting to check 
whether the current and charge densities are conserved in this case. The boundary conditions 
used in this work are 

iP{z)\l = iP{z)\r , -^^\l = ^—^\j^ (33) 

m{z) dz m(z) dz 

where L and R stand for left and right side of the heteroj unctions. This choice of the boundary 
conditions ensures that current and charge densities are conserved in the Hermitian system 
[30\ I40| . In this section we explore the possibility of this choice for the non Hermitian case 
discussed here. It is known from earlier works that charge density (say u) and current density 
(say j) for non Hermitian quantum systems with constant mass, obey the equation of continuity 
(for exact or unbroken VT symmetry) 



|^ + Vp-j = (34) 
only if the current and charge densities are redefined as 

00 = (f)*r](l) , j = i [^r](f) - (t>*V^^ (35) 

where is a linear, invertible, Hermitian operator, with respect to which the non Hermitian 
Hamiltonian H is pseudo Hermitian : 

ijt = n~^Hri (36) 

The interesting point to note here is that rj does not have a unique representation [17] . For VT 
symmetric potentials consisting of an even real part Vr, and an odd imaginary part Vj, viz., 

V{z) = VR{z)+iVi{z) or V{p) = VR{p)+iVi{p) (37) 

rj may be represented by the parity operator V [47\ I48j . For the non Hermitian VT Scarf II 
potential, viz., V{p) = —Vr sech^p — iVj sech p tanh p some definite forms of r/ are given in ref. 
With rj = V, and the forms given in ref. |49j . it is easy to check that for exact or unbroken 
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VT symmetry, the equation of continuity is obeyed in the transformed coordinate system p, 
where the problem reduces to one with constant mass, viz. eq. (j34p . Now, if one applies the 
inverse transformations of those in eq. (|1U|) . one can verify by straightforward algebra that eq. 
(j34p in the transformed coordinate system p (with constant mass) can be mapped to a similar 
equation of continuity in the original z-coordinate system (where the mass of the particle is 
dependent on its position) : 

§ + V.-J = (38) 

provided the charge density (w) and current density (j) in the original space are mapped to 
those {uj and j respectively) in the transformed space by 

Thus, if the charge and current densities {lo and j respectively) for non Hermitian VT symmetric 
quantum systems with position dependent effective mass are given a modified definition in 
accordance with eq. ()39p . then the boundary conditions in eq. (j33p ensure conservation of current 
so long as VT symmetry is unbroken or exact. However, with the spontaneous breakdown of 
this space-time symmetry, the current is no longer conserved. 



5 Conclusions and Discussions 

To conclude, in this work we studied a special form of semiconductor device consisting of a thin 
layer of VT symmetric material sandwiched between two normal semiconductors, such that the 
mass of the charge carrier (electron or hole) varies with position within the heterojunctions, but 
is constant outside. The mass m{z) and the real part of the potential Vr{z) are taken to be 
continuous throughout the device. We obtained the exact analytical solutions for the scattering 
states of a particle inside such a semiconductor device and also the refiection and transmission 
amplitudes, R and T respectively. Additionally, we also obtained explicit relations for current 
and charge densities for such a VT symmetric double heterojunction. 

The primary aim of this work was to extend the pdem formalism into the complex domain, to 
see if the spatial dependence of mass introduces any new feature in case of non Hermitian VT 
symmetric double heterojunctions. It is observed that at least two of the general features of VT 
symmetric potentials — viz., left-right asymmetry and blowing up of refiection and transmission 
coefficients at a spectral singularity — are preserved even for particles with pdem. The effect of 
the pdem is simply to introduce a non-linear component in the otherwise plane wave solution, 
within the heterojunctions. Another interesting feature we discussed here is the equation of 
continuity. For non Hermitian VT quantum systems with pdem, a modified definition of charge 
and current densities as given in eq. (j39p renders the conservation of current for exact VT 
symmetry. To the best of our knowledge, this is a new observation. 

Treating particles with position dependent effective mass because of varying doping concentra- 
tion in semiconductor devices of extremely small dimensions, is found to give better explanation 
to experimentally observed phenomena. On the other hand VT symmetric waveguides fabri- 
cated from iron-doped LiNbOa are also a reality [131 I50j . The time may not be far when a 
new generation of sophisticated, integrated devices are constructed based on the simple model 
discussed in this work. So our next attempt would be to see if the results of this study are 
typical of this particular example, or model- independent. 
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